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Abstract 
> 

oo 

O ■ particles (non-Abelian anyons) in an external magnetic field are addressed. We derive the 

(^ . A^-body Hamiltonian in the (anti-)holoniorphic gauge when the Hilbert space is projected 

onto the lowest Landau level of the magnetic field. In the presence of an additional 



The quantum mechanics and thermodynamics of SU(2) non-Abelian Chern-Simons 



harmonic potential, the A^-body spectrum depends linearly on the coupling (statistics) 
parameter. We calculate the second virial coefficient and find that in the strong magnetic 
Qj \ field limit it develops a step-wise behavior as a function of the statistics parameter, in 



contrast to the linear dependence in the case of Abelian anyons. For small enough values 
of the statistics parameter we relate the A^-body partition functions in the lowest Landau 
level to those of SU(2) bosons and find that the cluster (and virial) coefficients dependence 
on the statistics parameter cancels. 
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I. INTRODUCTION 

Identical particles with statistics continuously interpolating between Bose- Einstein and 
Fermi-Dirac statistics exist in two and one dimensions. Explicit implementations of these 
ideas can be found in the two-dimensional anyon model [|I| and in the one- dimensional 
Calogero- Sutherland models [0 where a different approach to statistics proposed by Hal- 
dane based on a generalized Pauli exclusion principle is realized. 

Contrary to the A^-body Calogero model, which is solvable, the A^-anyon spectrum 
is still unknown. However, a simplification arises when considering the anyon model in 
the background of an external magnetic field. By projecting this model onto the lowest 
Landau level, a procedure which is justified in the strong field-low temperature limit, a 
complete eigenstate basis, which continuously interpolates between the bosonic and the 
fermionic basis, can be found in the screening regime where the flux carried by the 
anyons is antiparallel to the external magnetic field B. More precisely, when the statistics 
parameter a = <f)/<f)o, which varies from a = to a = ±1, is such that a G [—1,0] if 
eB > 0, or equivalently a G [0, 1] if eB < 0. 

In this situation, the statistical mechanical properties of the anyon gas can be studied 
in a complete and explicit way 0] and they turn out to be quite similar to those of the 
Calogero model. In the thermodynamic limit both models are microscopical realizations 
of Haldane statistics @J^J§]. Recently, various conformal field theories have also been 
shown to implement exclusion statistics ^. 

It is well known that anyons can be thought of as bosons or fermions coupled to Abelian 
Chern-Simons gauge fields. An interesting generalization occurs when the gauge fields 
take values in a non-Abelian group and the particles carry internal degrees of freedom 
associated with a representation of this group. These models describing non-Abelian 
Chern-Simons particles have already been considered in several contexts. Verlinde |§] 
argued that they provide an explicit realization of non-Abelian braiding statistics, i.e. 
statistics corresponding to non-Abelian irreducible representations of the braid group M. 



Field theoretical implementations of such models were recently proposed [TQ] as Ginzburg- 
Landau Chern-Simons theories for Pfafiian (non Abelian) quantum Hall states. They 
generalize Abelian Chern-Simons field theories [|1^] for Laughlin (Abelian) quantum Hall 
states. 

Whether models of non-Abelian Chern-Simons particles restricted to the lowest Lan- 
dau level can be solved exactly (a natural possibility since they become effectively one- 
dimensional), and whether their thermodynamics yield a realization of exclusion statistics 
different from that corresponding to Abelian anyons, are questions of interest. In this pa- 
per, we address these questions by studying the simplest case of a non-Abelian symmetry. 



namely SU(2) Chern-Simons particles in the fundamental representation. The paper is 
organized as follows. 

In Section |ll| the Verlinde model [§ , a generalization of the Aharonov-Bohm Hamilto- 
nian in the (anti-)holomorphic gauge to the non-Abelian case, is introduced. We include 
contact 5'^{zi — Zj) interactions to enforce that the wave functions are regular when par- 
ticles approach each other, and, in addition, a confining harmonic potential to lift the 
lowest Landau level degeneracy. We redefine the wave function essentially as ip = Uip', 
where U (the non-Abelian generalization of Yli^j{zi — Zj)" for Abelian anyons with statis- 
tics parameter a) is a solution of the Knizhnik-Zamolodchikov equation. Our main result 
is the Hamiltonian acting on tp', which takes particularly simple form when ip' is (anti- 
)holomorphic, that is when the Hilbert space is restricted to the lowest Landau level of 
the external magnetic field. In the latter case we find a linear dependence of the A^-body 
energy spectrum on the statistics parameter, generalizing the known results for Abelian 
anyons. 

In Section |T| we revisit both the Abelian and non-Abelian anyon models to present 
another derivation of the same Hamiltonian acting on ip' as in Section ^ but now starting 
from free Pauli Hamiltonians. In this approach, contact interactions do not need to be 
introduced, instead the spin coupling plays a crucial role. We also show that in the 
redefinition of the wave function the external magnetic field and the magnetic fluxes 
carried by particles can be treated on the same footing. The latter point of view is more 
adapted for models that use mean field approximations [|12| . 



In Section |^ we discuss the 2-body problem for an arbitrary strength of the external 
magnetic field. The 2-body problem decomposes into two Abelian anyon problems, which 
allows one to calculate the second virial coefficient exactly. In the vanishing magnetic 
field limit our result reduces to that found previously by Hagen [^ in his comment on 



the paper by Lee |14]. In the strong magnetic field limit we find that the interpolation 
from Bose to Fermi statistics is different from the Abelian one. In the latter case 0, 
which yields Haldane's generalization of the Pauli exclusion principle, the second virial 
coefficient depends linearly on the statistics parameter. Contrary to that, we find that 
in the non-Abelian case the second virial coefficient develops a step-wise behavior as a 
function of the statistics parameter. 

In Section |V| we address the A^-body problem in the lowest Landau level. We show how 
to calculate the partition function taking into account properly degeneracies associated 
with internal isospin degrees of freedom. The A^-body energy spectrum for a given isospin 
is the A^-body SU(2) bosonic spectrum plus a a term linear in the statistics parameter. 
We relate the A^-body bosonic partition functions Z^j with total isospin I to the partition 
functions associated with Young diagrams. The latter partition functions, which turn out 



to be generalizations of the Schur functions, are introduced in subsection |V lj|. We propose 



systematic rules to calculate them. Collecting all these results, we can calculate the cluster 
and virial coefficients one by one. For small enough values of the statistics parameter we 
find, somewhat surprisingly, a cancellation of the dependence on the statistics parameter, 
and thus the same cluster and virial coefficients as for SU(2) bosons. The A^-body lowest 
Landau level thermodynamics of non-Abelian Chern-Simons particles in the entire interval 
of definition of the statistics parameter is yet an open question. 

We conclude in Section 0, and comment also on generalizations to other symmetry 
groups and relations to one-dimensional integrable models with inverse square interac- 
tions. In Appendix ^ some basic facts used in the paper on the irreducible representations 
of the symmetric group are collected. 

II. THE ABELIAN ANYON MODEL AND NON-ABELIAN VERLINDE 

MODEL 

A. Abelian case 

Let us begin with reviewing the formalism used in Q] to calculate the equation of 
state of Abelian anyons in the lowest Landau level of an external magnetic field. The 
dynamics of A^ anyons in the plane is described by a Aharonov-Bohm Hamiltonian in the 
background of a constant magnetic field B. Here, as in the sequel, the mass m will be set 
to 1 and complex coordinates notation Zi = Xi + it/i, Zi = Xi — iyi will be used. For our 
purposes, it will be convenient to work in the holomorphic (+) or anti-holomorphic (— ) 
gauge where the Aharonov-Bohm Hamiltonian takes the form 

^o^ = -E(V.^vf + vfvf) (1) 

i 

The covariant derivatives are defined as, 

\/f = ^,-^Kf^ Vf = a,-^ir| (2) 



with 



K = -^T.-^-^T.b-z^ K = o (3) 



K^.=^i:j^ + ^Y.bz, K=0 ' (4) 



The charge e of each anyon is coupled to the Aharonov-Bohm flux carried by the other 
anyons and to the external magnetic field B: therefore the couplings a = e0/27r and 
b = eB/2. 

We are interested in a feoson-based description of anyons, that is, the Hamiltonian (|T]) 
is acting on bosonic wave functions "^{zi, ..,zm]Zi, ..,Z]si). Then, the statistics parameter 
a G [—1, 1] is such that a = corresponds to bosons and a = ±1 to fermions. 

It is also possible, via a singular gauge transformation, to formulate the problem in 
terms of a "free" Hamiltonian, that is without anyonic Aharonov-Bohm fields, at the 
expense of introducing multivalued wave functions 

U+H+U-' = - Y.{dA + dA - 2bzA - h) (5) 



U-H^UZ^ = - Y.i^^^i + didi + 2hzid, + h) (6) 

i 

with U± satisfying the equations 

d,U+ = -J2 -^U+ d,U+ = (7) 

a 



diU- = J2 —U- d^U- = 



jVi ^i ^3 



which can be considered as the Abelian version of Knizhnik-Zamolodchikov equations [jT5 
with solutions 

[/+ = n(^. - ^.)-° V- = ui^^ - ^jT (9) 

i<j i<j 

Note that due to the multivaluedness of U, didiU ^ didiU . 
Hq have singular terms which arise from 

^^{h = d,{-) = 7f6\zi) (10) 

Zi Zi 

However, a potential accounting for repulsive contact interactions (which are introduced 
in order to implement the exclusion of the diagonal of the configuration space) and a har- 
monic well (which is introduced as a regulator to calculate thermodynamical properties) 
can be added to (0) 

V=J^XS\zi-z,) + Y.'^z,Zi . (11) 



If we choose A = 7r|Q;|, the contact term in V cancels exactly the singular terms in Hq 
{Hq) when a < (a > 0). The total Hamiltonian i^AB = Hq + V becomes 

^AB = - E (^^^^ + BA - 2b-zA " ^^.^i - b] + 2aY^ ^^ (12) 

HIb = -J2 i9^^^ + BA + 2bzA^ - ^% + &) - 2« ^ ^^ (13) 

The role of contact interactions in the context of anyons has been extensively discussed in 
the literature. Their relevance has been originally stressed in the study of soliton solutions 
in Chern-Simons matter systems |jl6| and anyonic wave functions in the background of 
an external magnetic field |T^ (see also [@). The issue was subsequently re-considered 



by several authors, in the perturbative treatment of the Aharonov-Bohm problem ||19|| , 
perturbative calculations of statistical and thermodynamical quantities [^, self adjoint 



extensions [^, etc. Notice that in the case under consideration the contact interactions 



are repulsive and that in addition, the orientation of the external magnetic field and the 
anyonic flux tubes are opposite, a case which do not support sohtons [^. On the other 
hand, this is precisely the case considered in Q since it allows for a physical meaningful 
lowest Landau level reduction. 

Then, if a; = 0, the ground state of H^^, which corresponds to 6 > and a < 0, is 
given by analytic wave functions while the ground state of H^^, associated to 6 < and 
a > 0, is given by anti-analytic wave functions. 

As we mentioned before, the harmonic attraction lifts the degeneracy of the ground 
state. In order to see this, it is convenient to re-define the wave functions as 

'^^{zi,..,zn;zi,..,Zn) = n^^P( ^-^ — ZiZi)ip^{zi, .., Zf^;zi, ..,Zm) (14) 

where uJt = V^^ + b'^- The Hamiltonian acting on i\)'^ is then 

if + = - ^ [diBi + BAi - {^t + b)'Zidi - {uJt - b)zAi - tOt) 

i 

H~ = -J2 {BiBi + BAi - i^t + b)ziBi - {uJt - b)zAi - c^t) 

-2«E(|^-^) (16) 



Acting on analytic and anti-analytic wave functions 

v^+ = n4' r = Yl^f (17) 

i i 

the Haniiltonians if ^ have the spectrum 

En = Nut + (E^^ - 2^^^ - ^))(^* - ^) ^^^) 



Em = Nut + (E^* + 9^(^ - 1))K + ^) (19) 

Notice that the exclusion of the diagonal of the configuration space, in view of (|^) and 
our discussion on the sign of a, is realized for the wave functions in the singular (s) gauge 

^^+ = n(^. - %)-" n ^ exp(-^z,^,) (20) 

i<j i 

r- = n(^. -^.)"n^"exp(-^z,^,) (21) 

i<j i 



which do vanish at coinciding points. 



B. Non-Abelian case: the Verlinde model 

Let us generalize the above construction to SU(2) non-Abelian Chern-Simons particles 
in the lowest Landau level of an external Abelian magnetic field. We are considering the 
case where A^ identical bosonic particles are in the isospin 1/2 /undamenta/ representation 
of SU{2). The wave functions of the system \Ef(2;i, .., z^, zi, . . . , zn) belongs to the tensor 
product Fi X ... X F^v where Fj is the two-dimensional space in which the SU(2) generators, 

T^='^ a = 1,2,3 2 = 1,.., AT (22) 

act. The non-Abelian generalization of the Aharonov-Bohm Hamiltonian in the holomor- 
phic and anti-holomorphic gauges is given by, 

^o^ = -E(V.^vf + vfvf) (23) 

i 

where we the covariant derivatives are defined as before but with non Abelian gauge fields, 

^:=^^E^^-^E^^^ ^i = (24) 



^^" = -*^E^?^ + *E^^^ ^.1 = ' (25) 

Here g is the Chern-Siinons coupling constant {g = 1/2ttk, in |]14[) and satisfies the topo- 
logical constraint g = 2/n with n an integer. 

Naturally, other choices such as Coulomb or axial gauges are possible. The connection 
among these different gauge choices is obtained via W~^HqW where W is not necessarily 
unitary. Unlike the Abelian case where W is known, for the non-Abelian case an explicit 
expression is available only in the 2-body case. 

The Hamiltonian acts on totally symmetric wave functions with bosonic interchange 
conditions, meaning that the wave function is symmetric under the interchange of both 
coordinates and isospin indices, 

y^ii,..,im,.-,«!v,«iv(^l5 ••' ^rn, ■■, Zl, .., Z]y) = '^ii,..,ii,..,i^,..,ij^{Zi, .., Zi, .., Z^, ••, Zn) (26) 

As in the Abelian case, there exists a Hamiltonian without anyonic gauge fields but acting 
on wave functions with non trivial boundary conditions. This Hamiltonian, U±Hq U^^, 
takes the same form than in the Abelian case where U± satisfies the non-Abelian Knizhnik- 
Zamolodchikov equations 

d,U+ = U+g Y, ^^ d,U+ = (27) 

r-parpa 

BiU^ = -U^g Y, _'_i diU- = (28) 

jjLi ^i - Zj 

As in the Abelian case, the Hamiltonians Hq have singular terms that can be elimi- 
nated by adding an appropriate potential V . Thus, as a non-Abelian generalization, we 
are led to consider if^g = H^ + V with 

V = \g\T,Y. T^T^S'iz. - z,) + Y ^-^^^^ (29) 



It can be shown that this potential corresponds to a repulsive interaction in the bosonic 
sector. We then obtain 

^AB = - E {9^^^ + OA ' 2bzA - "^z.z, - b] - 2gJ2 T^T^^^^ (30) 

HIb = -Y i9^^^ + BA + 2bzA^ ' ^z.z^ + b] + 2gY TtTf^^^ (31) 



and factorizing the Gaussian factor in the wavefunctions as in (p^) 

H^ = -Y, {^A + ^A - (^t + b)zidi - {uJt - b)zidi - ujt 

i 

-2^E^"^"(^^-^) (32) 

H' = -Y^ (didi + Bidi - {uJt + b)zidi - {uJt - b)zidi - t^t) 

i 

+2gYTtT^^-^) (33) 



Acting on analytic (anti-analytic) wave functions, H^ {H ) takes the form, 



where 



(34) 



(35) 



(li,M = gY.TtTf (36) 



It can be easily shown that the operator ^2 has eigenvalues, 

^N,i = I (/(/ + 1) - ^iv) (37) 

where A^ is as before the total number of particles and / is the total isospin. 
For analytic (anti-analytic) wave functions, the spectrum reads 

En,i = ATo^i + I E ^j + ^N,i I (^t - b) (38) 



N 

En J = Nut+\Y. ^J - ^N,i 1 {uJt + b) (39) 



III. THE ABELIAN AND NON-ABELIAN ANYON MODELS REVISITED 

In the standard presentation given above of the Abehan and non Abehan Verhnde 
anyon models, short range 6'^{zi — Zj) interactions have been added by hand to the non 
hermitian Hamiltonians expressed in the holomorphic or anti-holomorphic gauges. We 
now present another approach to derive the same Hamiltonians ( ^^3l) , which starts 
from free Pauli Hamiltonians, therefore implying spin coupling to the local magnetic field 
carried by the anyons. In addition, the eigenstates redefinitions with respect to the vortex 
(short distance) and the magnetic field (long distance) will be treated on an equal footing. 

In the singular gauge, let us start with the free A^-body Pauli Hamiltonians 

N 

HL, = -2Y.dA (40) 



TV 

Hi.ee = -'^T.9^^^ (41) 

i=l 

where the index ± refers here to the spin degree of freedom. 

If one considers the additional coupling to an external magnetic field, then, in the 
symmetric gauge, d —>■ d ~ bz/2 and d -^ d + bz/2. Of course one could choose other 
gauges for the B field, as the holomorphic or anti-holomorphic gauges, in which case in (^0]) 
d ^ d — bz,d ^ d, and in (0), d -^ d+bz, d -^ d, which yield basically the Hamiltonians 
d^J^) discussed above. However, we insist at this point on using the symmetric gauge, 
since it is the natural gauge to work with in the presence of the singular Aharonov-Bohm 
flux tubes. 

Indeed, the anyon model is defined via the non trivial monodromy of the A^-body 
eigenstates of -fffree 

IpireeiZl, 2^2, • • • , Zn] ^1, ^2, • • • , Zn) = e"*"^fc<' ^*='^(zi, ^2, • • • , Zn] Zi, ^2, • • • , ^A^) (42) 

where J2k<i ^ki is the sum of the relative angles between pairs of particles. As said before, 
ip{zi,Z2, . . . , Z]\[;zi,Z2, . . . ,zn) is by convention bosonic in the regular gauge with the 
statistics parameter a = for Bose statistics, and a = ±1 for Fermi statistics. 

Looking at (^) as a singular gauge transformation, one would obtain, in the sym- 
metric gauge, a A^-anyon Aharonov-Bohm Hamiltonian in the background of the external 
magnetic field with ^na J2i<j ^^{zi — Zj) interactions and TJ2ib shifts, induced by the spin 
up or spin down coupling to the local magnetic field of the vortices and the homogeneous 
background magnetic field. The parameter a represents as usual the Aharonov-Bohm flux 
carried by the anyons in units of the quantum of flux. 

10 



The short range (contact) interactions have to implement the exclusion of the diagonal 
of the configuration space (Pauli exclusion), and thus have to be repulsive. So, depending 
of the sign of a, the spin up Hamiltonian ( ^Ol) (a G [—1,0]) or spin down Hamiltonian 
(" e [0, 1]) have to be used. 

To materialize the short range repulsion in the eigenstates, one proceeds by redefining 



'i'{zi,Z2,...,ZN;Zi,Z2,...,ZN) = Y[ \^i " ^j T^^l^^l, 2^2, • • ■ , Zn] Zi, Z2, . . . , Zn) (43) 

i<j 

This is nothing but saying that the eigenstates do vanish as quickly as r^" when particles 
i and j come close together (again the =F sign has been chosen accordingly to the sign of 
a). At this point, (^J3^) together give 



V'free = Y[{Zk - Zi) "t/^h 
k<l 



a G 



-1,0] 



(44) 



V'free = Y[{Zk - ZlY'^). a G [0, 1] 



(45) 



k<l 



where both Y{k<i{.Zk — zi)'"" and Y{k<i{zk — zi)" are solutions of the holomorphic and anti- 
holomorphic Knizhnik-Zamolodchikov equations in the Abelian case. The non hermitian 
Hamiltonian acting on tp rewrites 



H^ 



N 
i=l L 



did-i - i-fziZi + -{zidi - Zidi) 






(46) 



H 



N 

-2E 



didi 



^ „ j ZiZi -r ^yZiOi ZiOi 






Zi Zj 



(47) 



The external magnetic field did not yet plaid any role in the eigenstates redefinition, 
thus the Hamiltonians ( ^ , ^7\ j expressed in the holomorphic or anti-holomorphic gauges 
with respect to the vortices, but in the symmetric gauge with respect to the external 
magnetic field. Also, there is no such singular gauge as (H2|) for a homogeneous mag- 
netic field. Let us however extract from the eigenstates the Landau exponential factor 
exp(±| J2i ZiZi). It should be considered on the same footing as ]!« \zi — ^jh" in (|43|) , as 
one can easily realize by considering the 2-dimensional identity 
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f dr'\n\f-7^\ = 7rr'^/2 (48) 

It means that a magnetic field can be regarded as the average of a distribution of vortices 



121. Let us redefine 



V'free = n(^fc ~ ^') " exp(--6 ^ ZiZj)?/'+ (49) 

k<l ^ i 



V'free = YliZk - Zl)" exp{-bJ2 ^i^i)^- (50) 

k<l ^ i 

where the prefactors z/'free = U'^ip± are now solutions of the holomorphic and anti- 
holomorphic Knizhnik-Zamolodchikov equations in presence of the external B field 

{d, + lz, + aJ2^—)U+ = {d, + lz,)U+ = (51) 



{d,-^z,-aY,-^)U-=0 id,-^Zi)U-=0 (52) 

In order for (^,^) to be physically meaningful, i.e. short distance vanishing eigen- 
states due to the repulsive vortices, and long distance exponential damping due to the 
magnetic field, one concludes that in (^UD, a < 0,6 > 0, whereas in (pDD, a > 0,6 < 0. 
In both cases, the local vortices carried by the anyons are antiparallel to the external 
magnetic field, i.e. a screening regime. The resulting Hamiltonian in the holomorphic 
and anti-holomorphic gauges acting on ip± narrows down to 

N _ _ F) — F) 

H+ = -2j2 [m - bzA] +2aJ2 ^^^ (53) 

N f) ft 

H- = -2Y: [B^d, + bzA] - 2a ^ ^^ (54) 

Any analytic function of the variables Zi is an eigenstate of the holomorphic Hamilto- 
nian (p3D, whereas a anti-analytic function of Zi is an eigenstate of (|5^). From ( ^pO] ), 
one infers the infinitely degenerate ground state with zero energyj] 



^The lowest Landau level spectrum has been shifted downward by the spin induced shift. 
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^: 



free 



n( 



Zi Z n 



-"l[4^expi--bJ2z,Zi), i,>0 



(55) 



^: 



free 



i<j 



Y[zl'exp{-bJ2ziZi 



> 



(56) 



where the orbital quantum numbers £i = 0, 1, ...oo are such that the eigenstates are 
symmetric. If one leaves aside the anyonic prefactors Y[i<j{zi — Zj)~" and Y{i<j{.Zi — Zj)"", 
the A^-anyon ground state is a symmetrised product of 1-body Landau ground states with 
orbital angular momentum ii, as for ideal one-dimensional bosons. 

To compute its thermodynamical properties, one has to regularize the system at long 
distance. Adding a harmonic well confinement J^i^^'^ZiZi to the Hamiltonians ( ^PD , 
an A^-body anyonic eigenstate in the lowest Landau level of an external magnetic field is 
still entirely characterized by a product 1-body eigenstates with a given orbital quantum 
number ii. 

li b > 0, b = +u!c where Uc is is half the cyclotron frequency, and a G [— 1,0]. If 
b < 0, b = —uJci and a G [0, 1]. In the presence of the harmonic well, the eigenstates are 
still given by (^,^), but now with cUc ^ i^t = 
rewrite as 



uol + cj2. It follows that (^^ should 



^free = n(^fc ~ ^') " exp(--U;i ^ ZiZi)^}^ 
k<l ^ i 



(57) 



V'free = n(^fc ~ ^^T ^^V^^Ty^t Y^ ZiZi)l\}. 
k<l ^ i 



{51 



to get the holomorphic and anti-holomorphic Hamiltonians in the presence of the harmonic 
well 



N 



H' 



-2E 

i=l 

+2«E 

i<j 



didi 



l^i + l^c _ o ^t-i^c 



-Zidi 



Zidi 



di - dj ut - uJc 



Z% Zj 



+ Y,{uJt- Wc) 



(59) 



TV 



H 



-2E 

-2a E 

i<j 



OiOi ;; Zidi Zidi 



di - dj LOt - LOc 



Zi Zj 



+ E(^< ~ tuc) 



(60) 
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These Hamiltonians are identical, up to a constant energy shii^, to the Hamiltonians 
( [T3|JTB| ) obtained above since, when acting on regular wavef unctions, didi = didi. 

The A^-body spectrum is nothing else but the sum of the 1-body spectra tl^ = {ut — 
^c) + ^i{ujt —ujc), shifted by the 2-body statistical energy —^N{N — l)(co't —uj^a. So, the 
virtue of the harmonic confinement has been to partially lift the degeneracy with respect 
to the £j's, and to dress the spectrum with an explicit a dependence 



Ejv = N{ujt - ^c) + (E^^ T i^N{N - l)a)K - ^c 



(61) 



Let us now turn to the non-Abelian case : The non-Abelian Knizhnik-Zamolodchikov 
equations become in place of (0) and ( |52D 










(62) 



(9, -^z,)f/- + f/- (7^:1^ = (^^-k)^" 







ji^i ^i ^i 



(63) 



One can proceed exactly in the same way as in the Abelian case, i.e. start from the free 
Pauli Hamiltonians H^^^^ in ( ^PD , eventually coupled to the external magnetic field, and 
redefine the free eigenstates according to 



i)iree{zi, Z2, . . . , Zm] ^1, Z2, ■ ■ ■ , Zn) = U^i'±{zi, Z2, ■ ■ ■ , Zn] Zi, Z2, ■ ■ ■ , Zn) 



(64) 



The non Hermitian Hamiltonian, in the presence of an harmonic well rewrites in the 
holomorphic and anti-holomorphic gauges as 



TV 



H^ 



-2E 



didi - 



-2^E^"^" 



ZiOi 

2 2 

dj - Bj _ Ut -^, 



Zidi 



+ E(^t - ^c) 



(65) 



N r 



H 



-2E 



didi — Zidi — Zidi 



+2gY^T^Tf 

i<j 



2 2 

di - dj LJt - LJc 



Z4 Zq 



EK 



ijJr 



(66) 



a generalization of (|59| , |60D which coincides with (|3^j33|) . Accordingly the spectrum of (|65| 
(respectively (|66|) ) acting on analytic (respectively anti-analytic) eigenstates is 
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En,i = N{ut - uj,) + (5] ^. ± n^j) {oot - ooc) (67) 

i 

where Q has been defined in (^7\). 

The cases 6 > 0, a G [—1, 0] and 6 < 0, a G [0, 1] are equivalent. So in the sequel one 
will assume without loss of generality 6 = c^c > 0, a G [—1, 0], i.e. the holomorphic gauge. 



IV. THE THERMODYNAMICS OF TWO NON-ABELIAN CHERN-SIMONS 
PARTICLES IN AN EXTERNAL MAGNETIC FIELD 

A. Virial expansion 

In this section, we study the thermodynamic quantities in which anyonic statistics 
manifests itself. For this purpose the virial coefficients a^, which result from the expansion 
of the pressure P in terms of the particle density p, can be used 

1 oo 

-\nE = pP = J2a,p' (68) 

^ k=i 

Using the cluster expansion, 

oo 

lnE = Y.hz'', (69) 

where z = e^^ is the fugacity, it is possible to express the /c-th virial coefficient in terms 
of partition functions Zj with j < k. Indeed, from 

oo 

E=Y: Z^z'^ (70) 

^=0 



it follows that. 



and 



hi = Zi , 

1)2 = Z2 — -Z^ , 

63 = ^3 - Z2Z1 + ^Zl , (71) 



ai = 1 , 

as = -263 + 46^ , (72) 
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where a^ = ak/V^~^ and hk = bk/b\. 

Considering a harmonic potential as a regulator, the thermodynamic quantities in a 
box of infinite volume V can be calculated from those in a harmonic well of frequency uj 
by using the following "thermodynamic limit prescription": in the thermodynamic limit, 
the cluster coefficients b^ are obtained from those in the harmonic well 6^ by (depending 
on the cluster coefficient order k) ||23| 

1 V 



{kp^uj'^yi^ A: 



d ' 



(73) 



where d is the spatial dimension. In two dimensions, in the presence of an external 
magnetic field, (|73[) becomes 

VpL , (74) 



kl3{ujt - ujc) 



where pi = ^ is the Landau level degeneracy per unit volume. 



B. 2-body Abelian anyon case 

Let us consider two non- Abelian Chern-Simons particles in an external magnetic field. 
This problem can be solved exactly since it decomposes into two 2-body Abelian anyon 
problems. We start by discussing the Abelian anyon problem paying special attention to 
the peculiarities of the strong magnetic field limit which are important for understanding 
the non- Abelian case. 

The spectrum for the relative motion of two Abelian anyons in an external magnetic 
field B and a harmonic well u is 

Enm = (2n + 1 + |m - a\)ujt - {m - a)ujc , (75) 

i^nm = e^("'-")'rl'"-°lL„(u;,zz) . (76) 

The spectrum and eigenstates are periodic with period 2 in a since m has to be chosen 
to be an even (odd) integer for boson (fermion) based anyons. So one can always restrict 
ae [-1,1]. 

In the thermodynamic limit cj — > 0, the spectrum and eigenstates rewrite as 

m> a : Enm = (2ra + 1)uJc (77) 

Ipnm = Z^-'^LniuJ.Zz) (78) 
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or 



m < a 



E„ 



(2n + 1 + 2(a - m))^;^ 



(79) 



V'ri 



^LJuJrZz) 



(80) 



The projection onto the lowest Landau level corresponds to n = 0. This implies that 
the wave functions are analytic (anti-analytic) for m > a {m < a). Notice nevertheless 
that the projection is not well defined at the bosonic end (when a -^ 0^ for boson based 
anyons or a — > — 1^ for fermion based anyons). Indeed, let us look at the ground state 
n = in the the boson based description. If a G [—1,0], the analytic ground state basis 

state belongs to this basis. However, if a G [0, 1], 

state is anti-analytic and 

0+, 



(78) is complete since then the m 

the analytic ground state basis is incomplete since the m 

has the energy which varies linearly with a, joining the ground state basis when a 

as represented in Fig. 1. 



^ 3 



m = 



_L 



_L 



_L 



-1 -0.5 0.5 1 

a 
FIG. 1. The lowest Landau level and the first excited state of the 2-body boson based anyon 

problem. 



The 2-anyon second virial coefficient in a magnetic field has been computed in ||T8[ {b 
and / stand for boson based and fermion based anyons, respectively) 
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b f ^T ( 1 

02 = — -F - tann x 



gX/-2xa_^N 1 



-a 



±^ + — (e^(H-«)_r 



sinh X cosh x 



^1) 



where x = Puc and At = y/2¥j3 is the thermal wavelength. It is depicted in the Fig. 2 
both for low magnetic fields (low x) and for high magnetic fields (large x). 

In the limit of the strong magnetic field, one obtains for boson based anyons 

1 



2pi 



-l-2a) ae[-l,0] 



(82) 



a^ 



2pL 



-l-2a + A(l-e 



— lax"- 



ae[0,l], 



^3) 



-0.5 



0.75 


1 


r^ 


1 


- 


0.5 


- 


r 




- 


0.?5 




/ 


6^_____ 


^^^ 





^^^^^^^^^ 


v^ 


^ a 


- 


0.25 


1 


V 


1 


_ 




a 



0.5 



FIG. 2. Second virial coefficient as a function of the statistics parameter a: (a) -p- for zero 
magnetic field, (b) ^ for x = 1, (c) ^ for x = 5, (d) ^ for x ^ oo. 

For fermion based anyons, the resulting expression is 

4 = -^ (l - 2a - 4e-'"(^+°)) , (84) 



As the exponential factor is only relevant for a ^ — 1, (|84D can be rewritten as 
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4 = -^ (l - 2a - 4e-2"(^+°)) , a E [-1,0] , (85) 

4 = ^{l-2a), ae[0,l], (86) 

Note that as far as x is finite, the virial coefficients are continuous functions of a. 
These results should be compared to the ones that would arise if one projects into the 
lowest Landau level at the very beginning, 

4 = -^(-l-2a) ae[-l,0], (87) 

4 = -— (-l-2a + 4) ae[0,l], (88) 

4 = -^(l-2a), (89) 

^Pl 

a2 has now a jump at a = while 02(^1) ¥" '^2(1)- These discontinuities are a direct 
consequence of the projection onto the lowest Landau level which ignores the 771 = 
(m = 1) state leaving the ground state -or in other words the vanishing of the gap- in 
the boson (fermion) based description for a = {a = —1). These states, when properly 
taken into account, smooth the discontinuities as can be seen in Fig. 2. 



C. 2-body non-Abelian problem 

Consider now the 2-body non-Abelian SU(2) case. The isospin is either 1 = 1 (triplet) 

Xi,i = l + + > Xi,o = ^(l + ->+!- + >) Xi,-i = |--> (90) 



or J = (singlet) 



Xo,o = ^(l + ->-!- + >)• (91) 



The wave function is 

'^l,mi = C(^l, ^2) Xl,mj , V'0,0 = (P^i^U Z2) Xo,0 • (92) 

As we have seen above, the action of the Hamiltonian on the direct sum of two sectors 
r7=i © r7=o is diagonal in each sector. The eigenvalues of the operator ^2 are given by 
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(|39|). Thus, with respect to the coordinate part of the wave functions, there is a "bosonic" 
sector for / = 1 and a "fermionic" sector for / = 0, with statistics parameters —(7/4 and 
35'/4, and degeneracies 3 and 1, respectively. It follows that the second virial coefficient 
can be expressed in terms of the second virial coefficients for Abelian anyons as 



, total 



34(-f) + 4(f; 



(93) 



This result is a generalization for B ^ to the one obtained in [|l^. Remember that 
< g < 2: in the fermionic sector, 3f7/4 can fall outside the interval of definition [—1, 1]. 
In this case, the Abelian second virial coefficient should be extended periodically, which 
means that (|93| ) can still be used but with 3(?/4 -^ Sg/A — 2. As a consequence, a cusp 
appears in the second virial coefficient for g = 4/3, as shown in Fig. 3. 

In the vanishing magnetic field limit, the second virial coefficient interpolates between 
its end values for bosons and fermions with two internal degrees of freedom when g varies 
from g = to g = 2. 



a2pL 



0.25 
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^^ A 
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\ 
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1 





-0.5 



-0.25 
4 



FIG. 3. Second virial coefficient for the non-Abelian problem as a function of the statistics 
parameter of the bosonic sector — | for (a) x = 10 and (b) x = 70. 
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In the strong magnetic field limit, the lowest Landau level second virial coefficient 
reads ii < g < 4/3, 

a2{9) = -T^ (94) 

and if 4/3 < g < 2, 

a,(^) = ^(l-2e-^K^-0) . (95) 

It is close to a step function, changing abruptly near g = 4/3 on an interval of width 
~ l/{j3ujc), from its bosonic end value — l/4pL to its fermionic end value l/4pL. 



V. THE THERMODYNAMICS FOR A iV-BODY NON-ABELIAN PROBLEM 

IN A STRONG MAGNETIC FIELD 

We now ask the question about a possible generalization of the 2-body analysis to the 
A^-body case. 

Let us ffist remind that in the Abelian case, from the lowest Landau level spectrum 
( PP|J55| ) one can deduce |Q an equation of state which coincides with Haldane's statistics 
thermodynamics for a gas of particles whose 1-body spectrum is reduced to a single level 
of energy E = ujc- The virial coefficients are 

an = (--)"-^-{(l + ay - a"} (96) 

Pl ^ 

The critical ffiling z/^ = — 1/a where the pressure diverges describes a non degenerate 
ground state with all the £j's null. In the singular gauge, 

N 

2^ 



i<j 



When a = —1, one recovers a Vandermonde determinant built from 1-body Landau 
eigenstates. Incidentally, in the Haldane statistics point of view, when a = —m, the 
non degenerate ground state coincides with the Laughlin eigenstates at the critical ffiling 
z/cr = l/m. 

A. Lowest Landau level non- Abelian case 

For particles with isospin 1/2, there is a one-to-one correspondence between the total 



isospin of the system and the Young diagram for the isospin wave function [|24[. The 
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corresponding Young diagrams may have at most two raws, of length y + / and ^ — I for 
total isospin / and will be d denoted as [(y + 1){^ — /)]. Since the total wave function 
is symmetric under interchange of particles, the symmetry properties of the isospin and 
coordinate wave functions are determined by the same Young diagram (see Appendix ^). 
The space of solutions of the A^-body problem is decomposed into a direct sum of 
sectors with given values of the total isospin: 
1 1 

-(g)...(g) - = r.iv ©r._iv„-| ©■••©r._iv_, ©■••©r,_iv \n]®---®Vj_n \n] (98) 

2 2 2 ^ ^ ~ 2 ^ '^ \ ~ \-~\ ~ \-~[ 

The number of times the subspace with total isospin / appears in the decomposition (|98[) 

is given by 

A^!(2/ + l) 
dNj - rf[(^+,)(^_,)] - ^|^j^^),(|_j), (99) 

which is also equal to the dimension of the representation [{\N + I){\N — I)]. 
The energy spectrum (|38| ) can be represented in the form 

En,i = {ut - uj,)nM,i{g) + <:r^(^i, . . . , ^at), (100) 

where 

TV 

E^°r^(£i, ...,^M) = Nujt + {ut - ^c) E ^j (101) 

is (up to a constant) the spectrum of A^ SU(2) bosons in a harmonic well of frequency 
{ojt — cUc) in the isospin / sector (s); the numbers £j are discussed in detail in the next 
subsection. This implies that A^-particle partition function is 

K 

2 



Zn= Y. e-^^^'-^''^^''''^^^ Z^^f"^ (102) 

where Z^°|°"*^ is the partition function of SU(2)-bosons determined by the spectrum ( |101| ). 



For (7 = 0, (|102|) reduces to 



f N 



Zn= E zly^- = Y.z%_,zl, (103) 

where Z^ is the usual A^-particle partition function for bosons without internal degrees 
of freedom. The latter equality in ( |103| ) follows from the fact that the grand partition 
function for isospin 1/2 bosons is the square of the grand partition function of bosons 
without internal degrees of freedom. 

The partition functions for SU(2) bosons Z^Y^^ ^^^ discussed in the next subsection, 
where these are related to partition functions associated with Young diagrams (see (|117|) ). 
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B. Partition functions associated with Young diagrams 

We start with a system of identical particles without internal degrees of freedom in a 
harmonic well. The single-particle energy spectrum is 

ei = {ut - uoc)^ (104) 

with a = 0,1,2,..., where we have left aside the constant ground state energy shift ujt 
which has no effect on the thermodynamics. 

The 2-body wave function of arbitrary permutation symmetry can be represented as 
a linear combination of symmetric and antisymmetric wave functions. This corresponds 
to decomposition of the set of all possible values of the quantum numbers (£i, £2), £±,£2 = 
0,1,2,... for two particles into 

ii < £2 and £1 < £2 (105) 

Correspondingly, the Boltzmann partition function is the sum of the partition functions 
for bosons and fermions {q = e~^('^'~'^=)) 



^Boltz 



qy 



where 



Zb 



Z^ + Zb 



q 



(106) 



(107) 



(l-g)(l-g2)' (i_^)(i_g2)- 

In the 3-body case, a wave function of an arbitrary permutation symmetry can be 
decomposed into a sum of wave functions with three types of symmetry, bosonic, fermionic 
and mixed. Note that there are two (identical) representations of the mixed symmetry and 
that these representations are two-dimensional[]. To discuss decomposition of all possible 
values of the quantum numbers {£1, £2, £3), for three particles, one should consider standard 
Young tableaux labeling all possible symmetries of wave functions (see Appendix ^. 

With each standard Young tableau we associate a £- Young tableau which is obtained 
by the change i ^ £i for all the numbers inside the standard Young tableau. In the 3-body 
case, we then obtain four £- Young tableaux 













£2 
£3 






£1 


£2 


£1 


^' 


£1 
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FIG. 4. £- Young tableaux for three particles 



■••The above decomposition of the wave function is a decomposition of the re^u/ar representation 
of 53 into irreducible representations. 
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The second and third tableaux are associated with the two standard Young tableaux for 
the states of mixed symmetry. 

All possible values of £i, £2, and £3 can be arranged as 

^1 < 4 < 4 , (108) 

^1 < 4 < 4 , (109) 

il<i2<i3, (110) 

^1 < 4 < 4 , (111) 



where (|108|) and (|111|) , corresponding to bosons and fermions, respectively lead to the 
partition functions 



1 



3\ 



(l-g)(l-g2)(i_g. 
^§ = (l-g)(l-g2)(l-g3) ' (112) 

Straightforward calculations show that ( |109|) and ( |110|) lead to the partition functions 
2g/[(l — g)(l — g^)(l — q^)] and 2g^/[(l — g)(l — g^)(l — g^)], respectively. We associate 
these contributions with the second and third £- Young tableaux, respectively, in Fig. ^. 
Thus the partition function of the Young diagram of mixed symmetry is 

2{q + q^) 

(i-?)(i-r)(i-r) 

The Boltzmann partition function rewrites as 

Zf°''' = ^ =Z^ + Zff + Z^ (114) 

as it should. 

We are now in position to propose general rules to write down partition functions 
associated with Young diagrams in the A^-body case. To this aim, we represent the 
partition function associated with a Young diagram Y^ with A^ boxes as 

Zy,=dy,^p^, (g)„=n(l-g'=). (115) 

[q)N k=i 

where dyi^ is the dimension of the irreducible representation of the symmetric group Sn 
corresponding to the Young diagram IV (equal to the number of standard Young tableaux 
for Yn)i and Pyj^ is a polynomial in q. Note that the partition functions for A^ Abelian 
bosons and fermions read Z\^ = l/[{q)j^] and Z}^ = q^ ''^ '^^ [{q) n]-, respectively. 
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To evaluate Pyj^iq), to each ^- Young tableau is associated a chain of inequalities for 
the integer numbers {ii < i'2 < • • • < ^at} in the following way: If the row containing 
the box with ij+i is located below the raw containing the box with ij, then ij < ij+i, 
otherwise, ij < C-j+i- Let si, S2, . . . be the positions of the signs '<' in this chain counted 
from the right. Then the contribution to the polynomial Pyj^ from a given £- Young 
tableau is q^^ ^^ . Finally, Pyj^ , associated with a given Young diagram, is the sum of the 
contributions from all £- Young tableaux for this diagram. 

For instance, for ( |109|) , Si = 1, and the contribution to Ppi] is q while for ( |11(]| ), Si = 2, 



and its contribution is q . For the 4 and 5-body cases, explicit expressions for Py^^ are 
given in Appendix 0. 

Notice a relation between the partition functions for Y^q and its conjugated Young 
diagram F^: if Py^ = YX=T CYr,{k)q'' then Py-^ = EIzT cy^(fc)g*^-"^''-^ 

As a consistency check, all possible partition functions in the A^-body case should sum 
to the A^-body Boltzmann partition function, generalizing ( |106| ) and ( |114| ) 

Using the partition functions above we checked using Mathematica, up to the 6-body case, 
that ( |116| ) is indeed valid. 

Now we consider SU(2) bosons. The decomposition of the space of solution of the 
A^-body problem into sectors with different isospins was discussed in the previous section. 
Taking into account the multiplicity 2J + 1 for isospin I states, the partition function for 
SU(2) bosons in the isospin / sector(s) can be related to the partition function associated 
with the Young diagram corresponding to isospin / as follows 

^bosons ^ ^21 + l)Z[(^^j)(^_^)j . (117) 

Note that Zr,N_^_J^.f^N_J^^ receives contributions from (ir(iv_^_jwiv_j-x, £- Young tableaux. This 
corresponds to taking into account all dr/N^J^,N_J■.^ sectors of given I in the decomposition 
of the space of solutions of A^ SU(2) bosons in calculating their partition functions. One 
can check, using ( |115| ) and the expressions for the polynomials given in this Section and 
in Appendix |B], that ( |117| ) leads to the correct expression for the total partition functions 
for N bosons ( [1031) . 



To conclude this section, we comment on the relation between the above partition 
functions and the Schur functions. Schur functions naturally arise for the system of 
non interacting particles whose 1-body spectrum is represented by a finite set of levels 
€1,62, ■ ■ ■ ,eM [ill- One defines xi = e"'^^^ . . ., xm = e"'^^". Then the Boltzmann 
partition function can be represented as a sum over all irreducible representations A of 
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^Af° ^(^1) • • • ) ^m) = (Xi + ■ ■ ■ + Xm) 



Y, d{X)Sx{xi,...,XM) , (118) 



A 

|A|=JV 



where S\{xi, . . . ,xm) are the Schur functions associated with the irreducible repre- 
sentations A of the symmetric group S^- The case discussed above corresponds to 
£k = {k — l){ujt — ujc) and M ^ oo. We thus observe that the partition functions dis- 
cussed in this section are a generahzation of the Schur functions to an infinite number of 
arguments corresponding to an infinite set of equally spaced 1-body energy levels. These 
functions can be used in studies of the statistical mechanics of systems with parastatistics 
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C. On the virial expansion for non-Abelian Chern-Simons particles in a strong 

magnetic field 

Unlike the Abelian case, an explicit expression for the A^-body partition functions is 
not available. Nevertheless, it is possible to construct iteratively the cluster coefficients 
b^ of order k for any k < N. The cluster coefficients in the thermodynamic limit are 
obtained as 

h = VkpL(3\im(ut-iuM (119) 

Using the partition functions (|102|) with ( |1171 ), we find the cluster coefficients of lowest 
orders 

5- = (-1 + ^2,0 + 3fi2,i)-^r^ r + - (120) 

&3 = K,0 + 3^2,1 - |(1^3 1 + ^3,|)) -^^ ^ + 

\ + ^2,0(1 - ^^2,0) + 6^2,1(1 - 7^2,1) - 

y z 4 

SI,. (1 - in,.) - 2n3,,(i - in.,)) j^^^ + ■ ■ ■ (121) 

In the thermodynamic limit, 6^ should behave as [iOt — ^cj^^ ■, i-e. for small uj bk should 
be proportional to the volume. It is easy indeed to verify that Q^j in (|37| ) is such that 
the correct thermodynamic limit is obtained. 

Moreover, Qjyj is such that, in the cluster expansion, the leading term (ut — ujc)^^ does 
not depend on g what we checked up to the fifth virial coefficient. This implies that the 
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cluster coefficients (and consequently the virial coefficients) do not depend on the statistics 
parameter. It would be interesting to understand a possible symmetry underlying this 
cancellation. 

The above results agree with the exact second virial coefficient analysis if g is suf- 
ficiently small, \g\ < 4/3. However, aX g = 4/3 there is a jump in the second virial 
coefficient. It follows that the A^-body spectrum (|38D should not be considered as the 
correct spectrum for arbitrary g. As we have seen in the 2-body case, a shift in the 
statistics parameter had to be done, which in turn lead to the jump in the second virial 
coefficient. We have not found a generalization of this procedure to the A^-body case. 

In analogy with the Abelian case, we can argue that the spectrum of anyons is indeed 
given by (|38|) if Vt^j are less than one. Under this assumption, it is easy to show that the 
cluster coefficients hk are independent of g when 

In this range the cluster coefficients are the same as those of SU(2) bosons. The behavior 
of the cluster (and virial) coefficients in the entire interval of the statistics parameter is 
yet an open question. 



VI. CONCLUDING REMARKS 

We have developed a formalism to study the thermodynamics of non- Abelian Chern- 
Simons particles in the lowest Landau level of an external magnetic field. We expect that 
the thermodynamics of non- Abelian Chern-Simons particles described by Hamiltonians 
of the type (|32|j33|) and (|65| , |66D with other groups, as well as other irreducible representa- 



tions, can be analyzed using the same ideas with minor modifications. The hope is that 
under some appropriate choices of the non- Abelian symmetry one can obtain non-trivial 
thermodynamics in the lowest Landau level. 

In this respect, an attractive possibility would be to use non- Abelian symmetries pro- 
posed recently in Ginzburg-Landau Chern-Simons theories for non- Abelian quantum Hall 
states |Ty] . An argument in favor of this choice is the new exclusion statistics thermody- 



namics found for non- Abelian quantum Hall states quasiparticles on the edge [B^ and for 



their bulk counterparts |2^. We plan to address this issue in future publications. 

Another observation is that the same form of the spectrum, namely a bosonic A^-body 
spectrum plus a linear term in the statistics parameter, arises in one- dimensional inte- 
grable models with inverse square interaction |^9|. In those models, the coupling to the 
Chern-Simons field is replaced by the coupling of the inverse square interaction. The 
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spectrum is then valid in the entire interval of definition of the coupling parameter. It 
follows that that the cluster and virial coefficients do not depend on the coupling param- 
eter in the entire interval of definition of the interaction parameter. It would certainly be 
rewarding to find a symmetry principle underlying this cancellation. 
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APPENDIX A: YOUNG DIAGRAMS AND IRREDUCIBLE 
REPRESENTATIONS 



For completeness we include in this appendix a brief description of the relation be- 
tween Young diagrams and representations of the permutation group of N particles Sn- 
In particular we show how wave functions of arbitrary permutation symmetry can be 
constructed using Young operators. 

The irreducible representations of S^ can be associated with the partitions of N in 
positive integers A^, 



Ai + A2 + ... + A^ 



N 



Ai > A2... > A, 



(Al) 



The different partitions (|A1|) can be depicted graphically by means of Young diagrams in 



which each number Aj is represented by a row of Aj cells. The partitions (Al) are usually 
denoted by [A] = [AiA2...Am] where the power is used to indicate the repeated appearance 
of the same integer. Thus for two particles we have, [A] = [2] and [A] = [1^] associated 
with the first two diagrams in Fig^, while for three particles we have [A] = [3], [A] = [21], 
and [A] = [1^] corresponding to the last three diagrams in Fig.^. 



FIG. 5. Young diagrams for two and three particles 
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The number of inequivalent irreducible representations of Sn is equal to the number 
of different Young diagrams. In our examples, 5*2 and S^ have two and three inequivalent 
representations respectively. Using Young diagrams it is also possible to determine the 
dimensions of such representations as well as to construct explicitly representations r''^^ (P) 
of the elements P G Sn and the basis vectors. To do that, a standard Young tablueax is 
defined in which the numbers from 1 to A^ are distributed among the boxes of a Young 
tableau in such a way that they increase from left to right along the same row and from 
top to bottom along the the same columns. Then the dimension of the representation 
is equal to the number of standard Young tableaux. For example, there is only one 
standard Young tableau for each of the diagrams in Fig 0, except for [21], and so they 
correspond to one-dimensional representations of the permutation group. The diagram 
[21] has two standard Young tableaux shown in Fig^, corresponding to a two-dimensional 
representation. 



1 


2 


1 


3 


3 




2 





FIG. 6. Standard Young tableaux for the [21] representation 

Since any element of Sn can be written as a the product of A^ — 1 transpositions Pi-i,i, it 
is sufficient to specify the matrices T^'^'{Pi-ii) representing these transpositions. The fol- 
lowing rules define the Young- Yamanouchi standard orthogonal representation [^ where 
superscripts r, t run over all standard Young tableaux for the Young diagram [A] . 

Diagonal elements 

1. Tl^J{Pi^ii) = 1 if in the tableau r, i and and i — 1 are on the same row. 

2. T}^{Pi^ii) = —1 if in the tableau r, i and and i — 1 are on the same column. 

3. T]^J{Pi-ii) = ±^, otherwise. Here d is the number of vertical and horizontal steps 
that one must take to move from i — ltoi, and the upper (lower) signs applies when 
the row containing the number i is above (below) the one containing i — 1. 

Non-diagonal elements 

1. Trf{Pi-ii) = (1 — ^)2, if the tableaux r and t differ only by a permutation of 
numbers i and i — 1 and where d is defined as above. 
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For instance, for the two-dimensional representation [21] we only need to find rt^^l(Pi2) 
and rPi](P23) (as P^ = P12P23P12 , P123 = P12P23 and P132 = Pfss)- The application of 
the rules above gives 

r[-i(p,,) = [ ^ M ' r[-i(p,3) = [ ^' 1^ ] • (A2) 

To construct the basis for the representation [A], A^! Young operators are introduced 



as 



u. 



[A] 
rt 



:§)^Erff(P)P. (A3) 



where fx is the dimension of the representation. These fx operators with a fixed second 
index transform into each other under permutations and form a basis for the irreducible 
representation F'"*'!. 

A basis more convenient in physical applications can be obtained by applying the 
Young operators to a product of "one-particle wave functions" in the following way. Con- 
sider a set of A^ orthogonal functions ipai'i), where i stands for the sets of variables on 
which the ipa depend, and form the product wave function \l'o, 

^o = M^)M^)---MN). (A4) 

One can then obtain a basis for the orthogonal representations by applying the Young 
operators to this function 

vj/ff = a;g'^o . (A5) 



[A] 
rt 



where P is understood as acting on the arguments of ^a. Then, the functions \1' 
with a fixed second index transform into each other under permutation and form a basis 
for a irreducible representation of the permutation group. The Young tableau r that 
corresponds to the first index, enumerates the different wave functions in a given basis 
and characterizes the symmetry under the interchange of arguments. The second index 
t enumerates the different basis for F^ and it can be shown that it characterizes the 
symmetry of the wave function under the interchange of the function ipa- 

One can apply the above general discussion to the case of the coordinate wave function 
of a 3-body system. Consider a coordinate wave function 0(123) = 0(ri, r2, r^) (eventually 
the product of single particle wave functions, 0(ri, r2, r^) = 0Q,(ri)0^(r2)0^(r3)). Then in 
addition to the totally symmetric and totally antisymmetric wave functions. 
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0[31 = ^(0(123) + 0(123) + 0(132) + 0(321) + 0(231) + 0(312) 
0[i'l = (^0(123) - 0(123) - 0(132) - 0(321) + 0(231) + 0(312)) , (A6) 

one obtains for the mixed symmetry wave functions for the two 2- dimensional represen- 
tation [21] 

05^1 = -i=(20(123) + 20(213) - 0(132) - 0(321) - 0(231) - 0(312)) , 

0ff' = ^(0(132) - 0(321) - 0(231) + 0(312)) , (A7) 

0g^ = ^(20(123) - 20(213) + 0(132) + 0(321) - 0(231) - 0(312)) , 
V 12 

0^2'^ = \ (0(132) - 0(321) + 0(231) - 0(312)) , (A8) 

For a system of bosons with two internal degrees of freedom, the total wave functions 
for a given total isospin / (symmetric under interchange of particles) can be written as 

^a/3 = E 0Sxi'/! , (A9) 

7 

where 0i^2 ^^*1 X7/3 ^-^e the coordinate and isospin wave functions of the permutation 
symmetry [A] = [{N/2 + I){N/2 - I)]. 

Again in the 3-body case, the eight-dimensional internal space can be written as the 
sum of a / = 3/2 four- dimensional space, associated with wave functions which are 
completely symmetric under the interchange of isospin indices, plus two 1=1/2 two- 
dimensional spaces associated with states of mixed symmetry. The basis of the / = 3/2 
subspace is therefore 

x^j = ! + ++>, 



2'2 



1/1 

X|,i = ^(1 + +- > +1 + -+ > +1 - ++ >) 



2'2 



X|,-i = ^(|--+>+|-+->+| + — >), 

X3,_3 = 1 > . (AlO) 

For 1 = 1/2 there are two two-dimensional representations, one generated by 

Xi i = ^^(2| + +- > -| + -+ > +1 - ++ >) , 



2'2 



v^ 



XiL = ^(2| - -+ > -I - +- > +1 + — >) (All) 



2' 2 
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and the other by 



xfl = ^i\ + -+>-\-++>) 



'2 

xl^=^(|-+->-| + — >)• (A12) 

The superscript in these formulas is related to the permutation symmetry. For a fixed 
total isospin projection, say Iz = 1/2, Xi i cind Xi i form a basis for a two-dimensional 

2'2 2'2 

[211 [211 

representation of S^. More precisely, they correspond to "^[a and '^2a\ with a = 1 or 

a = 2 having taken \E'o(l, 2, 3) = | H h >. Notice that since there are only two available 

states, + and — , a = 1 and a = 2 give the same results so that only two independent wave 
functions can be formed. The same argument applies to the Iz = —1/2 sector starting 
from^o(l,2,3) = | - +- >. 

According to (|A^), the total wave functions for the system of bosons are obtained by 
multiplying coordinate and isospin wave functions of appropriate symmetries. For the 
sector I = 3/2, one has four wave functions 

^3 = 03 „ fzi, Z2, Z3)X'i,mr ' (A13) 



,mj V'gjm/V^i? ^^' ^o/ A^,m/ ' 



where 03 ,„^(zi, Z2, z^) is a totally symmetric wave function (|A6|) . For two representations 
in the isospin 1 = 1/2 (mixed symmetry) sector, we obtain four basis wave functions 

i^^'^ =<P^''^X^P +0^.''^X?^ , (A14) 






2,ni7 ^21™/ 21™/ 2'™-f 2'™-f 

J t-xrir V^ y^ t~\ I 



(A15) 



where 0i J^ with a, 6 = 1, 2 are the same as 0„^ in (|A7| - |A8|) . 



APPENDIX B: POLYNOMIALS Pyjv(<5) FOR FOUR AND FIVE PARTICLES 



Following the rules of subsection [VB| , one obtains the following expressions for the 



polynomials PyS^) in Eq. (|115|) for four particles 

P^{q) = 1 

i^(g) = g + g2 ^ g3 

I^ = q' + q^ 

-.3 I „4 I ^5 



P^=q^ + q +q° 

^ = q' (Bl) 
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and for five particles 

^[4i](g) = g + g' + g' + g^ 
^[32] = g' + g^ + g' + g^ + g^ 

^[312] = g' + g' + 2g^ + q' + q' 
P[2H] =q^ + q'' + q'^ + q' + q' 
^[213] = g' + g' + g' + g' 
P[i5] = gio . (B2) 

Note that the sum of all the coefficients of -Py„(g) is equal to the dimension dy^ of the 
associated respresentation of Sn- 
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